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Symmetrization for fractional elliptic 
and parabolic equations and 
an isoperimetric application 


Yannick Sire* Juan Luis Vazquez t, 
and Bruno Volzone^ 


Abstract 

We develop further the theory of symmetrization of fractional Laplacian operators 
contained in recent works of two of the authors. The theory leads to optimal estimates 
in the form of concentration comparison inequalities for both elliptic and parabolic 
equations. In this paper we extend the theory for the so-called restricted fractional 
Laplacian defined on a bounded domain of BY with zero Dirichlet conditions outside 
of ft. As an application, we derive an original proof of the corresponding fractional 
Faber-Krahn inequality. We also provide a more classical variational proof of the 
inequality. 


1 Introduction 

In this paper we will develop further the theory of symmetrization for fractional Laplacian 
operators initiated in the papers [251 El EH], both in the elliptic and the parabolic set¬ 
ting, by extending it to a natural version of the fractional Laplacian defined on a bounded 
domain 0 of M. N which is known as the restricted fractional Laplacian. This research direc¬ 
tion combines classical themes in the study of nonlinear elliptic and parabolic equations, 

*Universite Aix-Marseillc, I2M, Centre de Mathematique et Informatique, Technopole de Chateau- 
Gombert, Marseille, France, E-mail: yannick.sire@univ-anm.fr 

^Departamento de Matematicas, Universidad Autonoma de Madrid, 28049 Madrid, Spain, E-mail: 
juanluis.vazquez@uam.es 

^Dipartimento di Ingegneria, Universita degli Studi di Napoli “Parthenope”, 80143 Italia. E-mail: 
bruno.volzone@uniparthenope.it 


1 



like symmetrization and accretive operators, with the recent interest in nonlocal versions 


of the diffusion operators, specially the fractional Laplacians. As an application of the 


obtained comparison results, we derive an original proof of the Faber-Krahn inequality 
(FKI) for such operators defined on the bounded domain LI. 

Before entering into the description of our results, we review in this introduction the 
necessary information about symmetrization, the elliptic-to-parabolic technique used to 
generate evolution semigroups, the precise definition of the fractional Laplacian operators 
and the relation among these topics. This constitutes a sort of review part of this paper. 

Definitions of fractional Laplacians on bounded domains. When working in the 
whole space domain W N there are several equivalent definitions of the fractional Laplacian 
operator (—A) 0 '/ 2 , 0 < a < 2, classical references being [39l 153]. The interest in these 
operators has a long history in Probability since the fractional Laplacian operators of the 
form (—A) CT / 2 are infinitesimal generators of stable Levy processes, see [HE 157]. Further 
motivation and references on the literature are given for instance in DfflEU- A particular 
definition that has been convenient for symmetrization purposes defines the operator for 
every given 0 < a < 2 as the trace of a suitable Dirichlet-Neumann problem via an 
extended potential function w that solves an elliptic equation in a upper half-space in 
T~i+ = x (0, oo) C M JV+1 . This is usually called Caffarelli- Silvestre extension [18]. It 
allows to reduce nonlocal problems involving (—A) CT//2 to suitable local problems (actually, 
a degenerate-singular elliptic equation), defined in one more space dimension. 

When we work on a bounded domain Li C M. N , things complicate because there are several 
options for defining the fractional Laplacian operator (—A) CT//2 . Two of them appear often 
in the recent literature. In our previous work m on symmetrization we have followed one 
of these approaches to define the fractional Laplacian as the Dirichlet-to-Neumann map, 
through an extended potential function defined in a cylinder C = Li x (0, oo) C R^ +1 , as 
was proposed in m and m ■ Zero values are assigned on the lateral boundary of C. We 
call this operator the spectral version of the fractional Laplacian on Li. Let us call this 
operator C\ (A stands for the Laplacian). This setting allowed us to derive in [65] the 
desired symmetrization results, which extend the standard symmetrization theory applied 
to elliptic and parabolic equations driven by the standard Laplace operator. But let us 
recall that there are remarkable restrictions on their validity in the form of conditions on 
the nonlinearities that are allowed in the equations. 

In this paper we will take the second usual approach to define (—A) 0- / 2 , which seems to 
be more natural in many applications. It consists in keeping the definition of fractional 
Laplacian in M. N but asking it to act on the null-extensions to M. N of functions u(x) defined 
in Li. So in principle we can use the most common formulation with a hyper-singular kernel 



( 1 . 1 ) 


on the condition that f(y) = 0 for y 0 D. Let us call this operator C. 2 - This option has 
been called the restricted Laplacian on a bounded domain | T9l [10] , but we will prefer the 
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name natural fractional Laplacian with Dirichlet conditions in this paper. The discussion 
on the relations and differences between the two types of operators on bounded domains 
is currently being investigated by several authors. Thus, Musina and Nazarov [44] use the 
name fractional Laplacian with Navier conditions for the spectral version, and fractional 
Laplacian with Dirichlet conditions for the restricted version. 

Here, we want to extend to operator C 2 the symmetrization theory we had developed for 
C\ in the papers [25] and [51]. This has an independent interest since there are subtle 
differences between the two operators, see mm- 

Symmetrization. Symmetrization is a very ancient geometrical idea that is used nowa¬ 
days as an efficient tool of obtaining a priori estimates for the solutions of different partial 
differential equations, notably those of elliptic and parabolic type. Since the topic is so 
well known let us only recall some facts that are relevant here. Symmetrization techniques 
appear in classical works like [Ml 07]. The application of Schwarz symmetrization to ob¬ 
taining a priori estimates for elliptic problems is already described in [66] and [12]. The 
standard elliptic result refers to the solutions of an equation of the form 

Lu = /, Lu = — ^2 di(dijdju ), 
hi 

posed in a bounded domain Ll C the coefficients { o - ij } are assumed to be bounded, 
measurable and satisfy the usual ellipticity condition; finally, we take zero Dirichlet bound¬ 
ary conditions on the boundary dLl. The classical analysis introduced by Talenti [541 [55i 
leads to pointwise comparison between the symmetrized version (more precisely the spher¬ 
ical decreasing rearrangement) of the actual solution of the problem u(x ) and the radially 
symmetric solution u(|x|) of some radially symmetric model problem which is posed in 
a ball with the same volume as D. Sharp a priori estimates for the solutions are then 
derived. Extensions of this method to more general problems or related equations have 
led to a copious literature. 

Elliptic approach to parabolic PROBLEMS. For parabolic problems this pointwise 
comparison fails and the appropriate concept is comparison of concentrations, cf. Bandle 
mm and Vazquez [58]. The latter considers the evolution problems of the form 

(1.2) dtu = AA(u), ri(0) = uo, 

where A a monotone increasing real function and uq is a suitably given initial datum 
which is assumed to be integrable. For simplicity the problem was posed for x £ M^, but 
bounded open sets can be used as spatial domains. The novel idea of the paper was to use 
the famous Crandall-Liggett Implicit Discretization theorem m to reduce the evolution 
problem to a sequence of nonlinear elliptic problems of the iterative form 

(1.3) -hAA(u(t k )) + u{t k ) = u(t k _i), k = 1,2, — , 
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where tk = kh, and h > 0 is the time step per iteration. Writing A(u) = v, the resulting 
chain of elliptic problems can be written in the common form 

(1.4) hLv + B(v) = f, B = A~ 1 . 

General theory of these equations, cf. [5], ensures that the solution map: 

T: f^u = B{v) 

is a contraction in some Banach space, which happens to be L 1 (fl). Note that the constant 
h > 0 is not essential, it can be put to 1 by scaling. In that context, the symmetrization 
result can be split into two results: 

(i) the first one applies to rearranged right-hand sides and solutions. It says that if two 
r.h.s. functions /i,/ 2 , are rearranged and satisfy a concentration comparison of the form 
fi ~< f 2 , then the same applies to the solutions, in the form B{v i) -< -B(w 2 )Q 

(ii) The second result aims at comparing the solution v of equation (11.41) with a non- 
rearranged function / with the solution v corresponding to /^, the radially decreasing 
rearrangement of /. We obtain that v is a rearranged function and B(v&) -< B(v), i. e., 
B(v) is less concentrated than B(v). 

This precise pair of comparison results can be combined to obtain similar results along the 
whole chain of iterations u(tk) of the evolution process, if discretized as indicated above. 
This allows in turn to conclude the symmetrization theorems (concentration comparison 
and comparison of L p norms) for the evolution problem (11.21) . This approach can be used 
in many different situations. In particular, it will be used below. 

Symmetrization for equations with fractional operators. The study of elliptic 
and parabolic equations involving nonlocal operators, usually of fractional type, is cur¬ 
rently the subject of great attention. Symmetrization techniques were first applied to 
PDEs involving fractional Laplacian operators in the paper [25] . where the linear elliptic 
case is studied: 

(1.5) (-A y' 2 v = /, 

The paper uses an interesting technique of Steiner symmetrization of the extended prob¬ 
lem, based on the Caffarelli-Silvestre extension for the definition of u-Laplacian operator. 
In [64] the last two authors of the present paper were able to improve on that progress and 
combine it with the parabolic ideas of [SS] to establish the relevant comparison theorems 
based on symmetrization for linear and nonlinear parabolic equations. To be specific, they 
dealt with equations of the form 

(1.6) d t u +{-A) a/2 A(u) = f, 0 < a < 2 . 

x For the definition of the order relation -<, see Section [3 
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Following the known theory for the standard Laplacian, the nonlinearity A is an increasing 
real function such that H(0) = 0, and we accept some extra regularity conditions as needed, 
like A smooth with A'(u) > 0 for all u > 0. The problem was posed in the whole space 
M. N . Special attention was paid to cases of the form A{u) = u m with m > 0; the equation 
is then called the Fractional Heat Equation (FHE) when rn = 1. the Fractional Porous 
Medium Equation (FPME) if m > 1, and the Fractional Fast Diffusion Equation (FFDE) 
if m < 1. Let us recall that the linear equation dtu + (— A) a ^ 2 u = 0 is a model of so-called 
anomalous diffusion, a much studied topic in physics. 

The results of [25] and [63] include a comparison of concentrations, in the form v# -< v, 
that parallels the result that holds in the standard Laplacian case; note however that no 
pointwise comparison is obtained, so the result looks a bit like the parabolic results of the 
standard theory mentioned above. Paper [64] considers both problems posed in the whole 
space and on a bounded domain. In the latter case the spectral fractional Laplacian is 
always chosen. 


2 Outline of results of the present paper 

We are interested in considering the application of such symmetrization techniques to lin¬ 
ear or nonlinear elliptic and parabolic equations with fractional Laplacian operators posed 
on a bounded domain, when the natural (i. e., restricted) version of fractional Laplacian 
is used. We denote the operator by £ 2 - 

Parabolic equations. To be specific, we want to treat evolution equations of the form 

(2.1) d t u + (-A) a/2 A(u) = /, 0 < ct < 2 . 

We want to consider as nonlinearity A an increasing real function such that H(0) = 0, and 
we may accept some other regularity conditions as needed, like A smooth with A'(u) > 0 
for all u > 0. The problem is posed in P, a bounded subset of M. N with smooth boundary. 
The parabolic result is developed in Section [5] and has to be compared with the results 
of papers [ 631 ) [65]. We will focus on the linear case A(u) = cu. This is the case that is 
needed in the isoperimetric application that we study in Section [6] 

Elliptic equations. The application of the method of implicit time discretization leads 
to the nonlinear equation of elliptic type 

(2.2) h(-/\f/ 2 v + B{y) = f 

posed again in the whole space P C M. N with zero Dirichlet boundary conditions; h > 0 is 
a non-essential constant, and the nonlinearity B is the inverse function to the monotone 
function A that appears in the parabolic equation (12.11) . The elliptic results are developed 
in Sections [3] and [4] and have to be compared with the results of papers [25l 163] where the 
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equation is posed either in or in O with operator C\. Note that the elliptic results we 
get cover the standard linear case where the term B(y) disappears and we set h = 1. 

A geometrical application. The Faber-Krahn INEQUALITY. As an application, we 
will use the symmetrization results to prove the Faber Krahn inequality for the fractional 
Laplacian operator on a bounded domain in both versions considered above. We recall 
that the FKI is a classical eigenvalue inequality, due separately to Faber m and Krahn 
[38] . based on a conjecture by Rayleigh in 1877, that can be stated as follows: 

Let Ll be a bounded domain in M. N and let B be the ball centered at the origin with 
V ol(Ll) = Vol(B). Let Ai(fi) be the first eigenvalue of the Laplacian operator, with zero 
Dirichlet boundary conditions. Then Ai(P) > Ai (B), with equality if and only if Ll = B 
almost everywhere. 

This is a classical result in the Calculus of Variations and proofs can be found in the clas¬ 
sical books like Chavel’s ED], see a recent proof in [15] . The question we want to address 
here is: will the result also hold for the usual versions of the fractional Laplacian operator 
(—A) 0 "/ 2 defined on bounded domains of M. N with zero Dirichlet boundary conditions? 

The answer is immediate in the case of the so-called spectral version of the Dirichlet 
fractional Laplacian , C\, since its eigenvalues, Afc(£i;D), are directly related to those of 
the standard Laplacian, Afc((—A);0), by the formula: 

(2.3) Afc(£ i; D) = (Afc((—A); D)) 0 '/ 2 . 

However, no simple relation like this one happens for the natural fractional Laplacian 
with the definition restricted type, L^. In Section [6] we will use our comparison results to 
present an original derivation of the fractional FKI. It does not make use of any variational 
interpretation, but only of some properties of the evolution process. The FKI can also be 
studied either by probabilistic or variational methods. For completeness, we also present 
a variational derivation, see more details in the mentioned section. 

Preliminary material and notation. In the paper we will use standard concepts and 
notation on symmetrization as fixed in [63! • We gather the main facts that we did not 
present here in the first appendix for the reader’s convenience. 


3 Elliptic Problem with lower-order term 


The case of the natural fractional Laplacian £2 will occupy our attention in this paper. 
We start our analysis by the following nonlocal elliptic problem with Dirichlet condition: 


(3.1) 


(—A) CT / 2 v + B(v) = / (x) in P, 
v = 0 in \ P, 
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where 17 is an open bounded set of R^, a € (0, 2) and / is an integrable function defined in 
17. We are interested in treating the case of a bounded domain 17 with the natural version 
of the fractional Laplacian. Exceptionally, 17 may be R^, but this case was treated in |64j . 
We assume that the nonlinearity is given by a function B : R + ->• R + which is smooth 
and monotone increasing with B( 0) = 0 and B'(y ) > 0. It is not essential to consider 
negative values for our main results, but the general theory can be done in that greater 
generality, just by assuming that B is extended to a function B : R_ —>• R_ by symmetry, 
B(—v) = —B{v). Note that we have changed a bit the notation with respect to equation 
(|2.2[) in the introduction, by eliminating the constant h > 0, but the change is inessential 
for the comparison results. 

For simplicity, all the restrictions of the fractional Laplacian operator will be denoted 
by (—A) CT / 2 . The underlying assumption is that such operator will be restricted to the 
ground domain of each boundary value problem where it is involved. 

The extension method. The Caffarelli-Silvestre method can be kept as extension to 
the whole 7~L + , with the following important proviso: the extension must act on the null- 
extensions to M. n of functions defined in 17. In view of this discussion, a solution to problem 
(|2.2|) is defined here as the trace of a properly defined Dirichlet-Neumann problem in the 
following way: 

-div^j, (y^Vw) = 0 

(3.2) M*,0) = 0 

1 dvo 

-lim y l ~ a — (x, y) + B(w(x, 0)) = / (x) 

K a o+ dy 

where 1~L + := R^ x (0,+oo), and K a is the constant K a := 2 1 " <T T(1 — |)/T(^), see [18], 
but such a value is not important here. 

3.1 Review of existence, uniqueness and main properties 

If problem (13.2|) is solved in an appropriate sense, then the trace of w over 17, Tro(u;) = 
w(-,0) =: v is said to be a solution to problem (13.11) . Note that the trace of w on the 
bottom hyperplane {y = 0} is the extension v of the function v defined in 17 by assigning 
the value zero outside of 17. This is what makes the difference with the case of the 
spectral Laplacian, where on the contrary the domain of the extended function w is the 
cylinder C = 17 x (0, oo) and w takes zero boundary conditions on the lateral boundary, 
T, = <917 x [0, oo), see [61]. In accordance with our choice of operator and in view of 
the iteration process that leads to the solution of the parabolic equations, we need only 
consider functions / that are restrictions to 17 of functions defined in the whole of R^. 

In order to make this more precise, we introduce the concept of weak solution to problem 


in R + , 

for x € R^ \ 17, 
for x G 17, 


7 



(13.211 . It is convenient to define the weighted energy space 


X a/2 (U + ) = Iw £ H} oc {U + ) '■ J + V 1 a \^x, y w(x,y)\ 2 dxdy < ooj , 
equipped with the norm 


(3.3) 


Mlx"/ 2 («+) := (yf H+ yl a \^ w (x,y)\ 2 dxdy^j 


\ 1/2 


For an open set E of R N we denote by H a ^ 2 (E) the classical fractional Sobolev space of 
order a/2 over E. We recall that for any u £ H a ^ 2 (M N ) there exists a unique a-harmonic 
extension w £ X rT / 2 (/H + ) of u to the half space 'H + , namely w solves 


(3.4) 


-div Xj j, ( y 1 a Vw) =0 in H + , 
w(x,0) = u(x) for x £ R N . 


Then we will write w 


Exty+{u). Moreover, for suitable functions u we have 

(_A ) a/2 u = lim y 1 -* *^-(x,y). 

bio- y —^ 0 + dy 


Now, in order to give a proper meaning of solution to a problem like (13.211 in a bounded 
domain fl, we define the space of all functions in X a ^ 2 (/H + ) whose traces over vanish 
outside of fl, namely 

(3.5) X^ 2 (U + ) = {w £ X a / 2 {U + ) : u4iv x{ o } = 0 in R N \ fi} . 

The domain of the natural fractional Laplacian (—A ) a / 2 is the space defined by 


(3.6) 


u(n) 


' H a / 2 (Q) 

< i?Jo /2 (fi) 

k K' 2 ^) 


if 0 < a < 1, 

if & = 1 

if 1 < a < 2, 


where H a i 2 (Q) and (O) are usual fractional Sobolev spaces, see [3D], and 
< / 2 (Si) = {» € ff 1 / 2 (Si) : jf 'EEix < 00} 
with d(x) = dist(x,dQ). It turns out that 

U(n) = {w\ nx{ 0 y.weX° / 2 ('H + )} 

(see ma for a detailed account on this question). Then we provide the following definition: 




Definition 3.1 Let 12 be an open bounded set o/I w and f £ -L 1 (f2). We say that w £ 
X^~(H + ) is a weak solution to (13.21) if TrQ(B(w)) =: B(w(x, 0)) € L 1 (0) and 

(3.7) / y l ~ a V x , y w-\7 X) yipdxdy + n a / B(w(x, 0)) <p{x, 0)dx = n a / f(x)ip(x,0)dx 

Jh+ Jn Jn 

for all the test functions ip € C 1 (fH + ) such that Tr R jv(<^) = 0 in \ 12. 

If w is a solution to the “extended problem” (13.21) . then the trace function v = Tr^(u;) 
will be called a weak solution to problem (13.11) . 

Remark 3.1 It is clear that if B(t) = ct for all t > 0 for some c > 0, then problem (13.21) 
becomes linear and the function v = Trn(w) belongs to the space 77(12). 

Concerning existence of solutions, their smoothness and L 1 contraction properties, we 
excerpt some known results from USES El! which can be extended for our more general 
nonlinearity B. For the regularity the reader may consult [IS IS USED EE]. 

Theorem 3.2 For any f £ L°°( 12) there exists a unique weak solution w £ X^- 2 (Cq) to 
problem (13.21) . such that Ttq(B(w)) £ L°°(12). Moreover, 

(i) Regularity: we have w £ C a (Cn) for every a < a if a < 1 (resp. w £ C 1,01 {Cq) for 
every a < cr — 1 if a > 1). Arguing as in [17] . higher regularity of w depends easily on 
higher regularity of f and B. 

(ii) L 1 contraction: if w, w are the solutions to corresponding to data f, f, the 

following L 1 contraction property holds: 

(3.8) [ [B(w{x,0)) - B(w{x,0))} + dx < [ [f(x) - f(x)\ + dx. 

Jn Jn 

In particular, we have that w > 0 in Cci whenever f > 0 on 12. Furthermore, if we put 
u := B(w(-,0)), then for all p £ [l,oo] we have 

IMIip(n) < H/llLp(f2)- 

(iii) For data f £ L 1 (12) the weak solution is obtained as the limit of the solutions of 
approximate problems with f n £ L 1 (12) n L°°( 12), f n —> / in L 1 , since then the sequence 
{B(w n (x,0))} n also converges in L 1 to some B(w(x, 0), and \\B(w(x, 0))||i < ||/||i, hence 
v n in uniformly bounded in L p for all small p. Property (ii) holds for such limit solutions. 


9 


Remark 3.3 (On some nonhomogeneous boundary value problems). Let e > 0. For our 
arguments it will be essential to consider problems with nonhomogeneous boundary values 
of the type 


(3.9) 


(—A)' 7 / 2 v + B(y) = / (x) in fl 


V = £ 


inR N \fl. 


In order to ensure the existence of a solution to this problem, we will associate to it the 
following nonhomogeneous extension problem 


-div X)l , ( y 1 a Xw) = 0 


in T-i + , 

for x € R N \ fl, 



—— lim y 1 a ^-(x,y) + B(w(x,0)) = f (x) 
0 + dy 


for x G fl. 


If / € L 1 (fl), we say that w € X rT ^ 2 ('H + ) is a weak solution to (13.101) if w — e G X^ 2 ('H + ) 
and w satisfies (13.71) . In such case, we say that v = Tr^rc is a weak solution to problem 
(13.91) . In particular, if B(t) = ct for all t > 0 and some c > 0, and v is the solution to the 
linear problem with homogeneous boundary data 

J (—A)' 7 / 2 v + cv = f(x) — ce in fl, 

[ 7=0 in R n \ fl, 

then v = v + e is the unique solution to the linear, nonhomogeneous problem ()3.9I) . 

Using the variational formulation to (13.101) . it is easy to prove that if e > 0, v £ is the 
unique weak solution to (13.91) and w £ is its extension solving (13.10|) , then 


w £ —^ w inL 1 ('H + ), 

Tt r nw £ Tr R jvi(;, inL 1 (M Ar ) 


where v, w solve m-m respectively. 

We warn the reader that the solutions of all the Dirichlet problems throughout the paper 
will be identified with their extension on the whole R N , whose values out of II will clearly 
depend on the boundary conditions considered. 

4 Concentration comparison for the extended problem 

Let us address the comparison issue. From now on, we will always assume that the right- 
hand side / is nonnegative. Our goal here is to compare the solution v to (13.ip with the 
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solution V to the problem 


(4.1) 


(-A) ff/2 V + B(V) = f# (x) in tl# 

V = 0 on r n \ n*. 


A reasonable way to do that is to compare the solution w to (13.2|) with the solution if to 
the problem 

— di v x ,y (y 1 _fr V^) =0 in % + 

2 ^ i ip(x, 0) = 0 for x £ 1^ \ tl# 

— — lim y 1 ~ a ^-(x, y) + B(ip(x, 0)) = f* (x) inti#, 

K a o+ dy 

where ip(x, 0) = V(x). According to [25], using the change of variables 2 = ( y/<j ) a , 
problems (13.21) and (14.21) become respectively 


(4.3) 


and 


(4.4) 


. d 2 


w 


—z‘ 


— A r w = 0 


dz 2 
w(x, 0) = 0 

( x > °) = (/ (x) - B(w(x, 0))) 


d 2 ib 

-z u ^Z - = 0 


dz 2 

V’ = o 

-Tj 7 (®,0) = cr u ~ 1 k cj (.f # (x) - B(t/;(x, 0))) 


in H + 

for x € \ tl# 

in 12, 

in Ti + 

for x <E R n \ tl# 
in tl#. 


where v := 2 (a — 1) fa. Then, the problem reduces to prove the concentration compari¬ 
son between the solutions w[x, z) and ip(x, z) to (|4.3D - (I4.4I) respectively. We now introduce 
the function 

(4.5) Z(s,z)= f (w*(t,z) —ij>*(T,z))dT . 

Jo 

Using standard symmetrization tools (see [25]) we get the differential inequality 


(4.6) 


- ( z v Z zz +p(s)Z ss ) < 0 
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for a.e ( s,z ) € (0,+ 00 ) x (0,+ 00 ), where p(s) = N 2 uj 2 ^ N s 2 2 / N . Obviously, we have 

(4.7) Z(0,y) = 0. 

A crucial point in our arguments below is played by the derivative of Z with respect to 
Due to the boundary conditions contained in (14.31) - (14.41) . we have for 0 < s < |fl 

(4.8) Z z (s, 0)>e a [ 8 (B(w*(T,0))-B(i/>*(T,0))dT 

Jo 

where 6 a := a a ~ 1 K a ] on the other hand, for s > |D| we have (dZ/ds)(s, 0) = 0. This is 
the novelty in the argument with respect to the spectral case treated before, where the 
boundary condition 

BZ 

— (|fl|,z)=0, Vz > 0 

made the real difference: the boundary conditions are imposed on y = 0 and are of two 
types. Contradiction must be obtained after taking into account the two possibilities. 


4.1 Comparison result for the elliptic problem 

We are going to obtain a comparison result for some linear and nonlinear B. Actually the 
nonlinearities considered here will allow to get a result which is weaker than the one for 
the linear problem, i. e., when B{t ) = ct for some c > 0, which is the only case that we 
are going to need in addressing the Faber-Krahn inequalities. We also point out that, in 
order to reach our goal, we will use a lifting-type argument of the symmetrized problem 

ED 

Theorem 4.1 Let v be the nonnegative solution of problem (ED posed, in a bounded 
domain with zero Dirichlet boundary condition, nonnegative data f £ L 1 (D) and B is 
smooth, concave, strictly increasing on M+ and such that B( 0) = 0. If V is the solution 
of the corresponding symmetrized problem ED. we have 

(4.9) v#{x)-<V(x). 

The same is true if LI = R N . 

Proof. In this case we pose the problem first in a bounded domain f l of R N with smooth 
boundary. We also assume that / is smooth, bounded and compactly supported in LI, 
since the comparison result for general data can be obtained later by approximation using 
the L 1 dependence of the map / *->• v. Let us choose e > 0, f e = / + B(e) and let us 
consider the solution V e to the problem 

f (~A) a/2 V E + B{V £ ) = ft {x) in LI# 

\ V e = e onR N \Ll#, 
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where ft = /# + B{e). By virtue of Remark 13.31 we have V £ = if e (x, 0) for all x € R N , 
where if £ is the solution to the problem 

- div Xi y (y 1-<T V'0e) =0 in Bt 

< ijj e (x,0) = e forx£R N \Q# 

—— lim y 1_cr ^(x, y) + B(ip £ (x, 0)) = / # (x) + B(s) in 12 # , 

Ka o+ dy 

which can be reduced to the problem 


(4.10) < 


Setting 


—2 


,d 2 lfe 


dz 2 
1 p £ {x, 0 ) = £ 
dif, 


- A x tp £ = 0 


Qz (x, 0) = (T a 1 n a (f* (x) + B(e) - B(i/j £ (x,0))) 

(s, z) = [ z ) - tp*(r, z))dr 
Jo 


in T~L + 

for xeR N \n* 
in S2 # . 


Z £ 


We will prove that 


(4.11) Z £ (s,z)< 0. 

To this aim, we first observe that Z £ satisfies inequality (|4.6|) . In particular, a big role is 
played by the property of the solution if £ : in particular, the level sets {x £ M. N : if £ (x, z) > t} 
are bounded because they are balls centered at the origin. Now we set 

Y e (s,z) = [ (B(w*(t, z)) — B(i/j*(t, z))) dr. 

Jo 

From (14.81) . we obtain 

BZ 

(4.12) —A{s, 0) > 0 ff Y e (s, 0) + 0 a B{e). 

dz 

By the strong maximum principle applied to equation (14.61) . , which is satisfied by a 
positive maximum of Z £ cannot be achieved at an interior point, hence it must be achieved 
either as s —> oo or at a boundary point (sq, 0) for some so > 0. The first option cannot 
hold since w*(s,z) —>• 0 while ip £ (s,z) —> £ as s —> oo. As for the second, we see that for 
s > | O | we have 

dZ 

-^-(s,0) = v*(s)-V e *( S ) = -£< 0 , 
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the function Z e (s,0) is strictly decreasing in [|fl|,oo), then it must happen that so £ 
(0, |0|]. Arguing as as in [64], using (14.121) . we can write 


(4.13) 


dz e ,_ „ r°,c. 


dZ F 


(s 0 ,0) > Q a J B'(v*(s)) : ^{s,0)ds 


= e n 


fS 0 

g(0)Z £ (s 0 , 0) + / [Z £ (s 0 , 0) - Z £ (s, 0 )]dg(s) 

■Jo 


> 0 


where g(s) := B'(v*(s)), which is impossible due to Hopf’s boundary maximum principle. 
Finally, by (14. IIP we have, for s > 0 


thus as e —> 0 


rs rs 

/ v*(r)dr < / V*{r)dT 

Jo Jo 


f v*{r)dT < f V*(r)dT. 

Jo Jo 


and the result follows. The case Q = has been solved in [6TI, Theorem 3.2], to which 
we address the interested reader. □ 


Remark 4.2 If we consider the linear case, i.e. when B(t ) = ct, for some c > 0, the 
proof of Theorem \J.1\ can be simplified, because Z £ = cY £ , moreover i/j £ = if + e where if 
solves (|4.4|h and V E = V + e. 


Remark 4.3 If B is nonlinear and satisfies the assumptions of Theorem \4-l\ the mass 
concentration comparison (14.91) is not enough to obtain the same result when passing to 
the evolution problem via Crandall-Liggett Theorem. 


A valuable property we are going to prove now is that equality in (14.91) implies that the 
domain f l of the initial problem (13.11) is actually a hall modulo translation. This kind 
of property is known to be essential for studying the equality case in the Faber-Krahn 
inequality, i.e., to prove that the ball is the unique minimizer of the principal eigenvalue 
of the classical Laplacian over all the sets of fixed Lebesgue measure (see for instance 
Kesavan [37]). 

Proposition 4.4 (Case of equality) Assume that we have an equality sign in (14.91) . in 
the sense that 

rs rs 

/ v*(a)dcr = / V*(a)da 

Jo Jo 

for all s £ [0, |0|]. Then Cl is a ball, i.e. Cl = Cl# (up to a translation of the origin). 
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Proof. Using the same notation as in Theorem 14.1 [ by the given assumption we have 
y(s,0) = 0 for all s € [0, |0|], and since the extensions are null, this equality holds for 
every s > 0. Then by Hopf’s maximum principle and (14.81) we find Y = 0, thus 

(4.14) w*(s,z) = tp*(s,z) for all s,z> 0. 


We divide the rest into some steps: (i) Here we argue as in [53] or in [32] • Recall that 
the function w(-,z) is smooth on M jV for any z > 0. Then let us fix z > 0, multiply both 
sides of the equation (|4.3D by the test function 


<Ph ( x ) = < 


1 if w(x, z) >t + h 

w(x, z) — t 


h 


0 


if t < w(x, z) < t + h 
if w(x, z) < t, 


and integrate over K jV . An integration by parts yields the identity 


\y x w\ z dx - z v ™ 

1 {xE.WL N :t<w(x,z)<t-\-h} J {xE.M. n : w(x,z) >t+h} 


L 


d 2 


dx 


—z 


7 

J{x 


d 2 w (w — t 


{x£S. N :t<w(x,z)<t+h} \ h 

Then if we let h —> 0 we find 


dx = 0. 


d_ 

dt 


f |V x ui| 2 dx = z v f 

J{x£.M. n : w(x,z)>t} J {: 


d 2 


W 


{x£M n : w(x,z)>t} 


dz 2 


dx 


thus using the second order derivation formula (see Appendix) we get 


-- / 

7 ^ J{x£R n : w(x,z)>t} 

Concerning the solution 7 to problem ( 14.41) . since it is spherically decreasing w.r.to x, the 
follow equality occurs 


|V x u;| 2 dx < z v 


riMv(t,z) g2 w * 


ds. 


-- f 

J{x£M. n : ip(x,z)>t} 


\V x if\ 2 dx = z l 


L 




dz 2 


ds. 


then using the fact that w*(s,z) = ijj*(s,z) (which implies fx w (-,z) = we have 


d_ 

dt 


lV a 


w 


dx < 


’ {ccER-^ :w(x,z)>t} 


-- 1 

J{x£M N : ip(x,z)>t} 


V x i/>\ 2 dx 
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Integrating between t and oo, we find 


/ \V x w\ 2 dx< / \V x ^\ 2 dx 

J {iGR w :w(x,z)>t} J {sGR^: ip(x,z)>t} 

then by the Polya-Szego inequality and (|4.14l) 

/ \S7 x vj#\ 2 dx < / \S7 x w\ 2 dx < / | S/ x i/j\ 2 dx 

J:w#(x,z)>t j- J{xEM N :w(x,z)>t} J{xEM N : ^(x,z)>t} 

(4.15) 

= / \W x w*\ 2 dx. 

J{x&L n :w #(x,z)>t } 

We conclude that for every t > 0 


|V x rc| 2 dx = 


iV^tc^l 2 dx , 


I {xGl* :w(x,z)>t} J -.w#(x,z)>t} 

which is the equality case in the Polya-Szego inequality. 

(ii) Now we notice that w(-,z) is analytic on the upper half space as a consequence of its 
representation formula. Indeed, the Poisson kernel for the extension operator L a is given 

by 


P(x, z ) = c 0 


,7V 


\ x \2 _|_ z 2^N+a 

(see e. g. EH), which is an analytic function on the upper half-space z > 0. This means 
that for every h{x) G L l (]R N ) the solution w(x,z ) of the elliptic equation with data 
w(x, 0) = h(x) will be analytic, since it is the convolution of h with P (with respect to 
the x variable). In fact, once we know that w G C°° in a certain subdomain, it will be 
analytic by classical results on solutions of elliptic equations with analytic coefficients (see 
for instance ESI, m, ESI). 

(iii) Moreover, each level set {x G M. N : w(x,z) > t} is bounded because w(-,z ) decays 
to zero as |x| —> oo. We may now use the equality case in the Polya-Szego inequality 
(see [14], [lSU]) to obtain that {x G R N : w(x,z) > i} = {x G : w#(x,z) > f} modulo 
a translation. Then all the level sets {x G M. N : w(x,z ) > t} are balls. The results also 
imply that for every fixed z > 0 the function w(-, z ) is radially symmetric up to translation. 

(iv) Finally, we take the limit z —> 0 and we conclude that u{x) is also radially symmetric 

(as a function defined in W N ). This means that the domain Q, which is the positivity set 
of u in M. N , must be a ball, i. e., II = 12^, up to a translation of the origin. □ 


The following result can be shown as in the proof of [6H Theorem 3.3] 

Theorem 4.5 (Comparison of concentrations for radial problems) Letv\,V 2 be two 
nonnegative solutions of problem posed in a ball Br{ 0), with R G (0,+oo] with zero 
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Dirichlet boundary conditions if R < +oo, nonnegative radially symmetric decreasing 
data /i, /2 G L 1 (i?fi(0)) and B(t ) = ct for some c > 0 and all t > 0. T/ien v\ and V 2 are 
rearranged, and 

(4.16) /i -< f <2 implies v\ -< V 2 ■ 


5 Symmetrization for the parabolic problem 


The theory of existence of weak solutions for the initial value problem 
(5.1) d t u + (-&) a/2 A(u) = f, 0 < cr < 2. 


with A{u) = cu m , and all c, rn > 0, has been addressed by the first author and collabora¬ 
tors in [40], and the main properties have been obtained. In particular, if we take initial 
data in L 1 n H~ s , then an //"^-contraction semigroup is generated, and the Crandall- 
Liggett discretization theorem applies. The construction and properties of the solutions of 
the evolution problem is thus reduced to an iterated application of the results obtained for 
the elliptic counterpart in the previous section. This has been carefully explained in [64] 
and is reviewed in Appendix III. Thus, using Theorem 17.41 below, we obtain the existence 
of a unique mild solution to the linear Cauchy-Dirichlet problem on the bounded domain 
P: 


(5.2) 


u t + (-a y/ 2 u = f 
u = 0 

u(x, 0 ) = uq(x) 


x g P, t > o 
x G M. N \ P , t > 0 
x G P. 


obtained as a limit of discrete approximate solutions by the ITD scheme. 

Concerning the application of symmetrization techniques to this type of parabolic prob¬ 
lems, we can employ Theorems 14.1114.51 and the arguments in the proof of [64 . Theorem 
5.3] in order to find the following result 


Theorem 5.1 (Concentration comparison) Let u be the nonnegative mild solution 
to problem (16.21) . 0 < cr < 2, posed in P, with initial data uq G L 1 (P), right-hand side 
/ G L 1 (P x (0, oo)). Assume thatuo G L l { P#) is rearranged, f(x,t ) G L 1 ( P# x (0,oo) is 
rearranged w.r. to any t > 0, such that 

U* <Uq 

and 

f#{;t)<J{;t) 
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for a.e. t > 0. Let v be the solution of the evolution problem 

vt + (—A ) a / 2 v = f(x,t ) x £ , t > 0, 

v = 0 x € l w 

x £fl* 


(5.3) 


x G M JV \ , t > 0 


v(x,0) = « 0 (x) 
Then, for all t > 0 we have 


(5.4) u#(\x\,t) -< u(|x|,t). 

In particular, we have ||it(-,£)|| p < ||«(-,i)|| p for every t > 0 and every p £ [l,oo]. 


Remark. The parabolic result only covers the linear equation, which is much below our 
original expectations, since the elliptic result covers indeed nonlinear equations. We do 
not know if the lack of nonlinear results is due to a failure of the expected form of the 
theorem, as copied from (15.11) . or is only due to a lack of technique. We remind the reader 
that in the case of the problem in the whole space we were able to prove the nonlinear 
symmetrization result when A is concave, and to show that the corresponding statement 
for A convex is false, so that the question was posed about which kind of statement could 
be true. Such question is broader in the present case. 


6 Application: an original proof of the fractional Faber- 
Krahn inequality 

Here we prove the validity of the Faber-Krahn inequality for the fractional Laplacian 
operator C 2 defined on bounded domains of M. N with zero Dirichlet boundary conditions. 
This operator appears often in theory and applications, and is known under the name 
of restricted fractional Laplacian, though we can call it the natural fractional Laplacian. 
Unlike the spectral Laplacian C\ , the spectral sequence {\k(C 2 '-> H)}fc is not directly related 
to the sequence of the standard Laplacian. However, it is known that the spectrum is 
discrete and given by a strictly increasing sequence {\ka/ 2 {ty = -^(£ 2 ;^)} (see e. g. 
m)- Our Theorem is then stated as follows: 

Theorem 6.1 We have 

(6.1) A W2 (P) > A W2 (P#) 

with equality if and only if LI = Ll#, up to translation. 

The proof we present here is completely elementary and uses neither the variational char¬ 
acterization (16.61) nor the nonlocal Polya-Szego inequality as in [13], but only the concen¬ 
tration comparison provided by Theorem 15.11 and the asymptotic definition of Ai ]0 y 2 (P), 
which can be derived by the decay rate of the parabolic problem ut + C 2 U = 0 . 
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Proof of Theorem \6.1\ Suppose that {4>k,a/ 2 ,ni x )} k are the (T 2 normalized) eigenfunc¬ 
tions of £2 and let us consider the function 


u(x, t) = e A Wa( n )* - 01 ^/ 2 ^ (x) 


solving the problem 


( 6 . 2 ) 


ut + (—A y/ 2 u = 0 
< u = 0 


k u(x,0) = 'lpl ) a/ 2 ,n( x ) 
By Theorem l5.ll we find, for all t > 0, 


x g n, t > 0 

x € \fl ,t > 0 

iSSl. 


«M)||2 < ||u(-,t)ll2, 


where v solves the problem 


(6.3) 


vt + (—A) a / 2 v = 0 
< v = 0 

. w (*.°) =^ta/2^ X ) 


Since we have 


x € , t > 0 , 

x € 1^ \ , f > 0 

x € n # . 


(6.4) 


| u (-, t )|| 2 = e- Xl ^ n)t 


and the expression for v can be given in terms of superposition, namely 

OO 

« 0 M) = X] C/ 2 ,! 1 >^, 0 -/ 2 , n#)L 2 ( n#)e" Afe ’ ,T/ 2 (n#)t V'fc, a/ 2 , «#(*) 


fc=i 


thus 


IK.0ll2 = 2l^W2 1 n^fcW2,n#)lia ( n# ) e 2Afc ’ CT/2(Q#)t 

k =1 

00 

< e- 2A w a ( n #) ^ |(< CTAn) 4, CT /2,n#)li 2 (n#) 

fc=i 

- P - 2 A 1 , CT/ 2 (n#)t|i ,# 112 _ P - 2 A lW 2 (n#)i 


This together with (16.4ji implies (16. ip . 

The case of equality. Analyzing the last list of inequalities that starts by ||?;(-,i)||| we 
conclude that in the case where A lj 0 y 2 (f 2 ) = ^i,o-/ 2 (^) we necessarily have 

= ll M (’> i )llL 2 (n) = e Ai,ct/2 ^° 
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so that we conclude that the coefficients of the Fourier expansion of vo = v(-,0) in terms of 
eigenfunctions are all zero but the first, in view of the known fact that the first eigenvalue 
\a/ 2 ^*) is simple. This means that vq{x) = 2 > q#, with a constant c > 0. By 

normalization we get 

( 6 - 5 ) i’ta/ 2,0( x ) = ^W2,n# 0 ) ■ 

But this is enough to apply the important Proposition 14.41 and obtain = 12^, and the 
result ends as before. □ 


6.1 A variational proof 


A direct proof of the FKI for our operator and similar can be based on the variational 
interpretation of the first eigenvalue, since it can be written as the minimizer of the 
Rayleigh quotient, where the local L 2 gradient energy norm is replaced by the Gagliardo 
seminorm (see Section [7] for some details). 

Suppose that {ipk,a/ 2 ,n( x )} k are the (£ 2 normalized) eigenfunctions of £ 2 - As already 
mentioned in the introduction, the proof of Theorem 16.11 is a direct consequence of the 
variational characterization of Ai j(T /2(f2). Indeed, we know by [48] that 


( 6 . 6 ) 


*W2 ( n ) = m in w 

ueH */ 2 {r n )\{o} 

u =0 on R n \Q 


/ / 

JR N JS. 1 


I u(x) - u(y) I s 
\x — y\ N+cr 


dxdy 


f 

Jn 


u 2 dx 


Then we could use the nonlocal (hence fractional) version of the Polya-Szego inequality 
(see for instance [45]) to see that replacing u with vft makes the Gagliardo seminorm in 
decrease, therefore (16.11) holds. Furthermore if equality occurs in (16.11) . the mini¬ 
mality of A li0 y 2 (f^) implies that ''/if , 2 ^(x) is an eigenfunction, but since the eigenvalue 
Ai,o -/2 (A^) is also simple, by normalization we have 


( x ) = V’l ,a/ 2 ,Q#( x ) 


which the same result in [45j shows to be possible only when fl = and V’i,<r/ 2 ,n = 

_LL 

^ 10/2 n U P i' 0 translation. 

More general version of the FKI. Actually, Brasco et al. m are able to establish 
a more general version of the FKI, which applies to a nonlinear variant of the fractional 
Laplacian, namely the fractional p— Laplacian. The main argument of the general varia¬ 
tional proof is the the use of a nonlocal Polya-Szego inequality, proved in [32], to estimate 
the first nonlinear eigenvalue. 
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Probabilistic approach. The Faber-Krahn inequality for the fractional Laplacian with 
Dirichlet data on a bounded domain of R N is stated with a hint of the proof based 
on probabilistic arguments as the last result, Theorem 5, of This means that the 
eigenvalues are also characterized in terms of the evolution, in their case the stochastic 
process of Levy type. Another proof with probabilistic flavor can be found in [9]. 


7 Appendices 

I. On symmetrization 

We gather here some basic information on symmetrization that can be useful to read the 
paper. We follow standard notations used in the literature, and we recall that we have 
presented a more detailed account in [63]. A measurable real function / defined on R N 
is called radially symmetric ( radial , for short) if there is a function / : [0, oo) —>■ R such 
that f(x) = f(\x\) for all x £ R N . We will often write f(x) = f(r), r = |x| > 0 for such 
functions by abuse of notation. We say that / is rearranged if it is radial, nonnegative 
and / is a right-continuous, non-increasing function of r > 0. A similar definition can be 
applied for real functions defined on a ball -Br(O) = {x € R^ : \x\ < R}. 

If P is an open set of R N and / is a real measurable function on P. We will denote by 
|-| the IV-dimensional Lebesgue measure. We define the distribution function p,f of / as 

Hf (k) = |{z <E P : |/ (x) | > k}\ , k > 0, 

and the decreasing rearrangement of / as 

f* (s) = sup {k > 0 : fif (k) > s} , s £ (0, |P|). 

We may also think of extending f* as the zero function in [|P|, oo) if P is bounded. From 
this definition it turns out that g,f* = /af (i. e., /, and f* are equi-distributed) and f* is 
exactly the generalized inverse of /i/. Furthermore, if ujn is the measure of the unit ball 
in R w and P^ is the ball of R N centered at the origin having the same Lebesgue measure 
as P, we define the function 


f* (x) = f*(u N \x\ N ) ,zep # , 

that will be called spherical decreasing rearrangement of /. From this definition it follows 
that / is rearranged if and only if / = /A 

Rearranged functions have a number of interesting properties. Here, we just recall the 
conservation of the L p norms (coming from the definition of rearrangements and the 
classical Cavalieri principle ): for all p £ [l,oo] 

II/IIlp(Q) = ||/*||LP(|o,n|) = \\f*\\LP(n#) j 


21 


as well as the classical Hardy-Littlewood inequality (see |34j ) 

r /* 1^1 r 

(7.1) / \f(x)g{x)\dx < / f*(s)g*(s)ds= / f#(x)g#(x)dx. 

Jn Jo Jn# 

• We will often deal with two-variable functions of the type 

(7.2) / : (x,y) £ Cq f (x,y) £ M 

defined on the cylinder := O x (0,+oo), and measurable with respect to x. Here H 
can be a bounded domain or M. N . For such functions, it will be convenient to define 
the so-called Steiner symmetrization of Cq with respect to the variable x, namely the set 
C# := x (0,+oo). Furthermore, we will denote by gf ( k,y ) and f* (s,y) the distribu¬ 
tion function and the decreasing rearrangements of (17.21) . with respect to x for y fixed, 
and we will also define the function 

f* (x,y) = f*{u N \x\ N ,y) 

which is called the Steiner symmetrization of /, with respect to the line x = 0. Clearly, 
/# is a spherically symmetric and decreasing function with respect to x, for any fixed y. 

• There are some interesting differentiation formulas which turn out to be very useful in 
our approach. Typically, they are used when one wants to get sharp estimates satisfied by 
the rearrangement u* of a solution u to a certain and it becomes crucial to differentiate 
with respect to the extra variable y (introduced in the extension process that is used in 
fractional operators) in the form 

vpO X,y)dx. 

u(pc,y)>u* (s,y)} Vy 

We recall here two formulas that have been already used in [25] and [64] 


Proposition 7.1 Suppose that f £ H 1 ( 0, T ; L 2 (H)) for some T > 0 and f is nonnegative. 
Then 

f*£H\0,T-L\ 0, |0|)) 

and if j {f(x,t) = f*(s,t)} \ = 0 for a.e. ( s,t ) £ (0, |fl|) x (0, T), the following differen¬ 
tiation formula holds: 


(7.3) 


'/ 0 , 2/)>/*0 


df 

,y) d V 


(x, y) dx 



dr. 


The second-order differentiation formula is 
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Proposition 7.2 Let f nonnegative and f £ W 2 ’°° (Cq). Then for almost every y £ 
(0, +oo) the following differentiation formula holds: 


[ d H 

+ ( 


% (*.»))' 


(x,y)dx = -tL j f* (r.y)dr — 

°y Jo Jf{x,y)=f*{s,y) |V*/| 




dTL N 1 (x) 
\ _1 

iN—l 


f(x,y)=f*(s,y) \^xf\ 


’f{x,y)=f*{s,y) \^xf\ 


dU N - v (x) 


• Mass concentration. We will provide estimates of the solutions of our elliptic and 
parabolic problems in terms of their integrals. For that purpose, the following definition, 
taken from [58], is remarkably useful. 


Definition 7.1 Let f,g£ L} (JRL N ) be two radially symmetric functions on . We say 
that f is less concentrated than g, and we write f -< g if for all R > 0 we get 



g(x)dx. 


The partial order relationship -< is called comparison of mass concentrations. Of course, 
this definition can be suitably adapted if /, g are radially symmetric and locally integrable 
functions on a ball Br. Besides, if / and g are locally integrable on a general open set Q, 
we say that / is less concentrated than g and we write again f < g simply if /# X g&, 
but this extended definition has no use if g is not rearranged. 

The comparison of mass concentrations enjoys a nice equivalent formulation if / and g 
are rearranged, whose proof we refer to [34] , m, IS9]: 


Lemma 7.1 Let f,g £ L 1 (fl) be two rearranged functions on a ball fl = Br( 0). Then 
f < g if and only if for every convex nondecreasing function : [0, oo) —> [0, oo) with 
<h(0) = 0 we have 

(7.4) f <h(/(x)) dx < f &(g(x))dx. 

J Q J 

This result still holds if R = oo and f,g£ Ll oc (K N ) with g —>■ 0 as \x\ —> oo. 


From this Lemma it easily follows that if / -< g and /, g are rearranged, then 

(7-5) II/IIlp(q) < llsllLP(n) Vp 6 [1, oo]. 

II. On Analyticity. In the analyticity argument of Proposition 14.41 we wanted to apply 
the results of [35]. Let us put 

N 

F(x, Z, U, Uj , Uj k) = Z U U( N+ 1)(7V+1) + U B 

1=1 
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(where subindexes indicate partial derivatives). Then, = z v for j = k = N + 1, 
= 1 for each j = 1,..., N and = 0 for j ^ k, thus 


v--\ dF 

(x, z, u , Uj,Ujk)CjCk = z^Cn+i + |Ci| 2 + ■ ■ ■ + I Cat | 2 > 0 

j,k =l ' 7 '' 

for all (x, z) G R^ -1-1 := $7, £ G R 77 " 1 " 1 \ (0,0). Then the equation 

z u w zz + A^u; = F(x, z, V x , z w, Vl jZ w) = 0 

is elliptic in D = R^ -1 " 1 . Since a solution w to such equation is C°° in My +1 and the 
function F(x, z,u,Uj ,Ujk) is analytic in ( z,Ujj ) G M_(_ x R Ar+1 we can apply the main 
Theorem in m and conclude that w(x,z ) is analytic in R^ - * -1 . 


III. On accretive operators and the semigroup approach 

Let X be a Banach space and A : D(A) C X —> X a nonlinear operator defined on a 
suitable subset of X. Let us consider the problem 


u'(t ) + A{u) = f, t > 0, 
u{ 0) = u 0 , 


where uo G X and / G L 1 (J; X) for some interval I of the real axis. For a wide class of 
operators, in particular the ones considered in this paper, a very efficient way to approach 
such problem is to use an implicit time discretization scheme that we describe next. 
Suppose to be specific that I = [0, T] (but this can be replaced by any interval [a, b] and 
the procedure is similar). The method consists in taking first a partition of the interval, 
say, tk = kh for k = 0,1,... n and h = T/n, and then solving the system of difference 
relations 


(7.2) 


Uh,k - Uh,k -1 

h 


+ A(uh,k) 


r(h) 

Jk 


for k = 0,1,... n, where we pose Uh,o = ^o- The data set j/^ : k = 1,..., n j is supposed 
to be a suitable discretization of the source term /, corresponding to the time discretization 
we choose. This process is called implicit time discretization scheme (ITD for short) of 
the equation u'(t ) + A(u) = f. It can be rephrased in the form 

^ h,k — Jh(Ah,k —1 T h/j, ) , 

where the operator J\ = (/ + A-4) -1 , A > 0 is called the resolvent operator , / being 
the identity operator. Therefore, the application of the method needs the operator A to 
have a well-defined family of resolvents with good properties. When the ITD is solved, 
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we construct a discrete approximate solution {uh t k} k - By piecing together the values Uh,k 
we form a piecewise constant function, Uh(t ), typically defined through 

(7.3) u h (t) = u hik if t € [(k - 1 )h, kh\ 

(or some other interpolation rule, like linear interpolation). Then the main question 
consists in verifying if such function Uh converges somehow as h —>• 0 to a solution u 
(which we hope to be a classical, strong, weak, or other type of solution) to problem (17.11) . 
To this regard, we first choose a suitable discretization | fj, h> | in time of the source term 
/, such that the piecewise constant interpolation of this sequence produces a function 
f( h \t) (defined by means of (17.31) 1 verifies the property 

II — /IIl 1 (o,t ; x) “^0 as h —>• 0. 

By means of these discrete approximate solutions we introduce the following notion of 
mild solution: 

Definition 7.2 We say that u € C((0,T);X) is a mild solution to (17.11) if it is obtained 
as uniform limit of the approximate solutions Uh, as h ^ 0. The initial data are taken in 
the sense that u(t) is continuous in t = 0 and u(t) uq as t —> 0. Besides, we say that 
u £ C((0, oo);X) is a mild solution to (17.11) in [0, oo) if u is a mild solution to the same 
problem in any compact subinterval I C [0, oo). 

In order to state a positive existence result, we need to restrict the class of operators 
according to the following definitions. 

Definition 7.3 Let A : D(A) C X —> X be a nonlinear, possibly unbounded operator. 
Let R\(A) be the range of I + XA, a subset of X. 

(i) The operator A is said accretive if for all A > 0 the map I + A.4 is one-to-one onto 
R\(A) C X, and the resolvent operator J\ : R\(A) —>■ X is a (non-strict) contraction in 
the X-norm (i. e., a Lipschitz map with Lipschitz norm 1). 

(ii) We say that A satisfies the rank condition if R\{A) D D(A) for all A > 0. In 
particular, the rank condition is satisfied if R\(A) = X for all A > 0; in this case, if A 
is accretive, we say that A is m-accretive. 

We are now ready to state the desired semigroup generation result, that generalizes the 
classical result of Hille-Yosida (valid in Hilbert spaces and for linear ^4) and the variant by 
Lumer and Phillips (valid in Banach spaces, still for linear *4), and provides the existence 
and uniqueness of mild solutions for problems of the type (17.11) in the case / = 0: 
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Theorem 7.3 (Crandall-Liggett) Suppose that A is an accretive operator satisfying 
the rank condition. Then for all data uo £ D(A) the limit 

(7.4) S t {A)u 0 = lim (J t/n (A)) n u 0 . 

n—too ' 

exists uniformly with respect to t, on compact subset of [0, oo), and u(t) = St(A)uo £ 
C([0, oo) : X). Moreover, the family of operators {Sf(.4.)} t>0 is a strongly continuous 
semigroup of contractions on D(A) C X. 

Using a popular notation in the linear framework, we could write St(A)uo = e~ t,A uo, and 
because of this analogy formula (17.41) is called the Crandall-Liggett exponential formula 
for the nonlinear semigroup generated by —A. The problem with this very general and 
useful result is that the X-valued function u(t) = St(A)uo solves the equation only in 
a mild sense, that is not necessarily a strong solution or a weak solution. Though it is 
known that strong solutions are automatically mild, the correspondence between mild and 
weak solutions is not always clear. For the FPME this issue has been discussed in detail 
in [MEZ]. 

In addition, the Crandall-Liggett Theorem result can be extended when we consider 
nontrivial source term /, according to the following result 

Theorem 7.4 Suppose that A is m-accretive. If f £ L 1 ( 0,oo;X) and uq £ D(A). Then 
the abstract problem m has a unique mild solution u, obtained as limit of the discrete 
approximate solution Uh by ITD scheme described above, as h —>■ 0; 

u(t) : = lim u h {t ), 
h—> 0 

and the limit is uniform in compact subsets of [0, oo). Moreover, u £ C([0, oo);X) and 
for any couple of solutions u\, U 2 corresponding to source terms f±, f 2 we have 

||ni(t) - u 2 (t )||.y < | \ui{s) - u 2 (s)||x + ||/i(t) - f 2 (r)\\ x dT 

for all 0 < s < t. 

There is a wide literature on these topics, starting with the seminal paper by Crandall and 
Liggett [24], see also [23] and the general reference [5]. These notes are based on Chapter 10 
of the book [60], cf. the references therein. The last formula we have mentioned introduces 
the correct concept of uniqueness for the constructed class of solutions. Characterizing 
the uniqueness of different concepts of solution is a difficult topic already discussed (with 
positive results) by Benilan in his thesis |6]. 
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8 Comments and extensions 


• We have only proved results on parabolic comparison based on symmetrization for the 
linear case. The elliptic results can be applied to nonlinear equations but still have severe 
restrictions. It will be interesting to know how much is true for nonlinear functions B and 
A in the respective equations. This question was partially addressed and solved for the 
spectral fractional Laplacian in [63], and the limitations to the generality of the results 
were also shown to be necessary, the symmetrization result was false for concave B or 
convex A of power type. 

More generally, we would like to know if there is an approach that ensures comparison 
results of some symmetrization type valid for quite general nonlinearities, as it happens 
in the non-fradtional case, cf. [59] . 


The variable coefficient case. As a future direction, we are interested in the following 


= 0 in R N \Q 

%w\ y=0 = J 


A(x ) 0 

0 1 

Here the matrix A{x) is supposed to be IT 1, °°(M Ar ) and uniformly elliptic with lower 
constant A > 0. 


problem 

r div( y l ~ 2s 

(8.1) 

{ 


[ -y 

where 

B(x) 


It is a by now well-known fact that the spectral powers of div(A(x)V), i.e. ^div(A(x )v)) ; 


for s G (0,1) in a bounded domain H can be described as the Dirichlet-to-Neumann 
operator of a suitable extension in a cylinder C = H x M + (see for instance m for a 
detailed account). The previous problem (18.11) is a variant of this extension but in the 
whole M N . The Dirichlet-to-Neumann operator in this case is not explicitly identified; 
however, we believe it is a natural possible extension of the problem we considered in 
this paper. The idea here is to develop the techniques produced in the present paper to 
handle variable coefficients, having in mind an isoperimetric inequality. Indeed, a FKI 
was proven in terms of the first eigenvalue by means of 


Ai(fl) > AAi(D # ). 


The aim here is to prove such a result for the following problem: let C s be the Dirichlet- 
to-Neumann operator associated to the problem (18.11) defined on D. It is obvious that 
C s has discrete spectrum {Afc^jfTp It is not clear how to use a variational approach to 
deal with this operator since it does not seem obvious that this operator is associated 
to a norm in WL N satisfying a Polya-Szego inequality. However, the parabolic approach 
developed in the present paper seems promising. 
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